Electron release of rare gas atom clusters under an intense laser pulse 
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Calculating the energy absorption of atomic clusters as a function of the laser pulse length T we 
find a maximum for a critical T*. We show that T* can be linked to an optimal cluster radius 
R* . The existence of this radius can be attributed to the enhanced ionization mechanism origi- 
nally discovered for diatomic molecules. Our findings indicate that enhanced ionization should be 
operative for a wide class of rare gas clusters. From a simple Coulomb explosion ansatz, we derive 
an analytical expression relating the maximum energy release to a suitably scaled expansion time 
which can be expressed with the pulse length T*. 



After a basic understanding of the mechanisms gov- 
erning atoms and molecules subjected to an intense laser 
pulse |], ||, analogous studies on clusters pioneered 
by Rhodes H and Ditmire have appeared over the 
last years with a recent spectacular culmination in the 
demonstration of deuterium fusion in clusters ||. Most 
of these studies do focus on the situation after the laser 
pulse, namely on the abundance and kinetic energy spec- 
tra of electrons and ions. Some discussion has been de- 
voted to the question if the expansion of the cluster is 
driven by hydrodynamics or by a Coulomb explosion. 
Only very little attention has been paid to this type of 
dynamics in the time domain ^|. This is even more 
surprising since the time scales involved show that the 
expansion of the nuclei occurs on the same time scale 
as the pulse lengths which can be chosen, namely some 
10 to 1000 fs, or roughly 10~ 3 atomic units (which we 
will use hereafter). Apart from the nuclear motion and 
the pulse length T energy absorption from a laser pulse 
and subsequent ionization and fragmentation of the clus- 
ter involve two additional time scales, the optical cycle 
2ir/u> = 0.055a.u. for the typically used Titan-Sapphire 
laser of 800 nm wavelength, and the period of the bound 
electrons, which is of the order (hydrogen) of 1 a.u.. 
We will work with peak intensities between 10 14 — 10 16 
W/cm 2 . 

In the following we will demonstrate that the seemingly 
complicated process of energy absorption and fragmen- 
tation in the laser pulse can be split into three different 
phases, an 'atomic ' phase I, a 'molecular' phase II, and a 
relaxation phase III. Phase I lasts for a time To after the 
pulse has begun and is characterized by boiling off elec- 
trons through multiphoton or tunneling ionization, hence 
we have termed it 'atomic' phase. We define it to last un- 
til every second atom in the cluster has lost one electron, 
or equivalently until the probability of loosing an elec- 
tron in an atom has reached p = 1/2. This probability 
is calculated from a Krainov tunneling rate |g] where, 
however, the instant electric field is formed by the laser 
and eventually already existing charged particles in the 
cluster. 



Up to To we may assume that the atoms/ions have not 
moved yet. The second, molecular phase is characterized 
by Coulomb explosion of the cluster. During this phase, 
as we will show below, the cluster expands to a critical 
radius R* which optimizes the energy absorption. Phase 
III finally, until the end of the laser pulse and beyond, 
sees a relaxation of the system and the full fragmenta- 
tion of the cluster proceeding. The existence of these 
phases follows from a careful analysis of our numerical 
results. The relevance of the phases is underlined by the 
time which is spent under phase II. This time turns out 
to be instrumental for relating the electron release quan- 
titatively to the laser and cluster properties, as will be 
shown below. 

To simulate the process of energy absorption numeri- 
cally we have developed a quasiclassical model for small 
rare gas clusters. The nuclei are treated completely clas- 
sically, with the initial configuration defined by minimiz- 
ing the pairwise Lennard- Jones interactions Elec- 
trons bound to an atom or ion at position R are charac- 
terized by an effective binding energy 

E b = E^ tom + V to UR) (1) 

where the exact atomic binding energy E^ tora is shifted 
by Vtotab the sum of the potentials from the laser field 
and all other charged particles except the atom/ion the 
electron is bound to. Ionization from such a bound state 
is accomplished via tunneling along the direction f of the 
instant force at position R , r = V V(R)/\VV(R)\. The 
time-dependent tunneling action along f reads 

1(f) = f 2 V2(V t otai(r) + VAtom(r) - Ebj dr (2) 

with the classical turning points r% determined by 
VtotaiM + V^tom (r) - E b = 0. From I(t) we get the 
tunneling probability P(t) — exp(— 27(f)) and finally the 
tunneling rate 

Ht) = 7^P(t) (3) 
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with the classical Kepler period Tk of an orbit with bind- 
ing energy E^ tova . For each time step dt, a random num- 
ber z is compared to the probability w(t) dt for ionization 
during this time step. If w dt > z, the electron is born as 
a classical particle and placed at the outer turning point 
Ti obeying total energy conservation. From then on, this 
electron follows Newton's equations, and the next bound 
electron can be ionized. Hence, strictly sequential ion- 
ization is enforced. 

The interaction between two particles with charge Q\ 
and Q2 and position vectors f\ and r*2, respectively, is 
described with a smoothed Coulomb potential 

V _ QlQ2 ,.s 

V Softcore — r— — — , 

V (r-i - r 2 ) z + ax{Qx) + a 2 (Q 2 ) 

where the a$ are charge-dependent soft-core parameters. 
For electrons we used a(— 1) = 0.1, while the ionic m are 
chosen such that the potential minima for each ion always 
coincide with the quantum mechanical binding energy. 
This choice prevents artificial classical autoionization. 

The model allows us to follow the full time-dependent 
evolution of the cluster with all interactions for a long 
time (10 5 a.u.) to investigate the influence of the clus- 
ter expansion during the laser pulse on its energy ab- 
sorption. Although it implies, e.g. for xenon clusters, 
to propagate up to 200 charged particles, the computa- 
tion can be handled with moderate resources due to the 
crucial simplification which arises from treating bound 
electronic motion not explicitly. 

After the pulse is over the electron release from the 
cluster is a typical observable which changes as a func- 
tion of the pulse length T as shown in Fig. |l| for a Xei6 
cluster in comparison with the corresponding electron re- 
lease (i.e., ionization) of a Xe atom. The energy content 
of the laser pulse £ = F(t) 2 dt has been kept constant 
which means that the peak intensity Fq of the pulse with 
amplitude F(t) = Fo sin 2 (nt/T) cosuit decreases with in- 
creasing pulse length T according to F§ <x 1/T. As a 
reference for this energy normalization we chose a pulse 
with F = 0.16 a.u. and a pulse length of 20 optical cy- 
cles. One sees that the ionization of an atom increases 
towards shorter pulse lengths T or equivalently, higher 
peak intensity. Indeed, atomic ionization depends on the 
peak intensity Fq rather than on the pulse length T which 
is obvious if the electron yield is dominated by sequential 
ionization depending exponentially on Fq via the Krainov 
rate but only linearly on the pulse length. The os- 
cillations in the single atom case are due to the atomic 
shell structure. 

For the cluster the situation is quite different: for short 
T (high intensity Fj 2 ) the electron release/atom follows 
that of an isolated atom. In fact, it is even slightly lower 
than in the atomic case, due to the much larger Coulomb 
field of a multiply ionized cluster which has not signifi- 
cantly expanded. However, the electron release increases 
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FIG. 1: Average atomic charge of Xeie (•) an d Xe (o) as a 
function of pulse length. The lines are to guide the eye. 

again and reaches a maximum for some optimum pulse 
length T* at considerably smaller peak intensity. Hence, 
the cluster expansion plays an important role for the en- 
ergy absorption, in contrast to the atom for which this 
degree of freedom does not exist. A comparison of time 
scales shows, that indeed the Coulomb explosion of the 
ions in the cluster happens on the same time scale as the 
pulse duration (10 3 — 10 4 a.u. or equivalently some 10 
to 100 fs). Hence, the dependence of the electron release 
on T points to the spatial expansion of the cluster which 
may in turn exhibit a maximum electron release for a 
certain cluster radius R. We define R in terms of the 
averaged distance between two ions in the cluster, 

/ N \" 2 

m=^^fif{\Ri-Rj\ 2 }) ■ (5) 

First we assess the influence of the size of the cluster on 
the electron release under the reference pulse of 20 field 
cycles. The size of the cluster is varied preserving its 
shape by scaling the ionic positions Rf = Ai?° with a 
factor A compared to the equilibrium positions R®. As 
can be seen from Fig. a critical value of the mean 
interionic distance, R* = X*Rq exists, where the ioniza- 
tion yield shows a maximum. The position of R* hardly 
changes upon variation of the laser frequency. The ioniza- 
tion yield, however, increases with increasing frequency: 
this is due to the smaller ponderomotive oscillation am- 
plitude at higher frequencies, which leads to increased in- 
teraction between quasi free electrons driven by the laser 
field and those still well in reach of the cluster ions. 

The mechanism responsible for the existence of R* 
was first identified in diatomic molecules under the name 
CREI or enhanced ionization (ENIO) jt0[ [0} with the 
(linear) laser polarization parallel to the molecular axis. 
It might seem astonishing that the cluster also exhibits 
ENIO although there is no preferred axis which could 
align with the polarization axis. This is even more sur- 
prising since no enhancement was found for diatomic 
molecules if the polarization is perpendicular to the 
molecular axis or if the laser is circularly polarized. How- 
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FIG. 2: Electron release of Xei6 for fixed nuclei. Results are 
shown for for a pulse length of T = 55 fs and the frequencies 
w = 0.055 a.u. (•), uj = 0.075 a.u. (*), and uj = 0.11 a.u. (□). 
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FIG. 3: Electron release of Xei6 with fixed nuclei and circular 
polarization. 



ever, a distinct feature of ENIO is the insensitivity to 
changes in the laser frequency which we also find in the 
cluster (Fig. |^). This fact, together with the relation 
of T* to the critical radius R* as presented below pro- 
vide sufficient evidence that intense laser field dynamics 
of clusters is structured by ENIO as is the corresponding 
dynamics of molecules. For clusters, ENIO is even more 
general since there is no restriction with respect to the 
polarization of the laser: The direction of the axis for lin- 
ear polarization does not matter and ENIO also occurs 
for circular polarization as shown in Fig. [|. 

Our findings exemplified here for Xei6 have been con- 
firmed by extensive calculations for a number of clusters 
of 8 to 30 atoms for the elements Ne, Ar, Kr, and Xe. 
These calculations clearly demonstrate that ENIO plays 
an important role for small rare gas clusters under in- 
tense laser fields with quantitative consequences as we 
will see next. 

One could think that the relation of T* and R* is di- 
rectly given by the Coulomb explosion mechanism. The 
latter links indeed R* to a certain time interval r, but 
not to the entire pulse length T*. The reason lies in 
the existence of the different phases as described in the 
beginning and schematically shown in Fig. 

Only phase II drives the relevant Coulomb explosion, 
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FIG. 4: Sketch of phases I, II and III during the pulse (see 
text) 



namely the onset of the cluster expansion. For the maxi- 
mum electron release this time interval ends if the critical 
radius R* is reached at half the pulse length T*/2 when 
the laser pulse intensity has its maximum. The interval 
begins, however, only at To when each atom in the clus- 
ter has on average a 50% probability of being ionized. At 
this time the cluster still has its equilibrium radius Rq. 
Hence, we get r = T*/2 — T * as the relevant time during 
which the cluster expands from Rq to R*. 

Having identified phase II and its time interval r as the 
one which controls the maximum number Q* of released 
electrons, we can use the dynamics of Coulomb explosion 
to derive a quantitative relation between Q* and the op- 
timum pulse length T*. To this end we assume that the 
ionic motion in r can be described by an averaged ionic 
charge which is proportional to the averaged final charge 
of each atom in the cluster, i.e. 



Qi = aQ*/N, 



(6) 



where N is the number of atoms in the cluster. Further- 
more we assume that in phase II the main kinetic energy 
goes into the expansion of the cluster without changing 
its shape. Using the same parameterization as before, 
we write now for the time-dependent position Ri(t) of 
an atom or ion Riit) = A(t)R®, where R® is the equilib- 
rium position of the atom before expansion, i.e., at time 
t = T . The kinetic energy K = M/2J2 l (dR i (t)/dt) 2 
reads then 



K = 



AI 



N 



To- 



(7) 



lo has form and units of a moment of inertia and rep- 
resents the influence of the shape of the cluster on its 
kinetic energy during the expansion. The potential en- 
ergy V = Y.i >3 QiQj/\Ri{t) - R](t)\ simplifies to 



iV 



V = \- 1 (aQ*/N) 2 £ = a 2 V a /X. (8) 



i>j=l 



The differential equation for the expansion in terms of 
X(t) is obtained via the energy balance Kit) + V(t) — 
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FIG. 5: Expansion time (numerical data) as a function of 
(Ko/Vo) 1 ^ 2 /(A) and linear fits (see text). Two different en- 
ergy normalizations were used: Fq = 0.16 a.u. (solid line) and 
Fo — 0.25 a.u. (dashed line), both at a frequency of ui = 0.055 
a.u. and a pulse length of T — 55 fs. •: Ari6, o: Ar2o, 
Ar 2 5, O: Ar 30 , *: Kri 6 and x: Xei 6 



E = V(Tq), where at time To before the expansion the 
kinetic energy of the atoms is zero. With the help of 
Eq. (uv and Eq. M it can be written in the form 



f = a[(l-A-> o/ 7o] 1/2 . 



(9) 



Eq. H can be solved analytically by separation of variables 
to yield 



t(A) - T 



K Q 
V a 2 



1/2 



/(A) 



(10) 



where we have set A = 1 for t = T) and /(A) := 
l) + ln(\/A — T+VA). For the maximum energy 
release the critical radius R* — \*Rq should be reached 
after time i(A cr it) = T*/2. This is the desired relation 
between the static ENIO mechanism at R* and its dy- 
namical effect during the cluster expansion in the time 
interval r = T*/2 - T *. 

The proportionality factor in Eq. determines the 
fraction of the final charge by which the expansion from 
Ro to R* during phase II is effectively driven. If phase 
II is indeed the crucial dynamical time span which uni- 
versally controls the electron release we expect a to be 
the same for all types of clusters we consider, indepen- 
dent of the atomic element or cluster size. Under this 
assumption, we predict from Eq. (|ic| ) a linear relation 
between the expansion time r and {Kq/Vq) 1 ^ 2 /(A). In 
Fig. U we can see that this prediction is actually very 



good: shown are the expansion times r as a function of 
the cluster-dependent values of (Kq/Vq) 1 ^ 2 /(A) for dif- 
ferent clusters. A linear fit to the data yields a = 0.38 
and a = 0.37 for energy normalized pulses at F = 0.16 
and T = 0.25, respectively. The correlation coefficient is 
in both cases higher than 0.99. Hence, a is the same for 
different clusters, and it is almost the same for different 
energy normalizations of the laser pulse. 

To summarize, we have shown that the enhanced ion- 
ization mechanism is operative for small rare gas clus- 
ters over a wide range of parameters. Moreover, from 
a careful analysis of the Coulomb explosion process, we 
conclude that energy absorption and subsequent ioniza- 
tion of the cluster proceeds in a very similar way for dif- 
ferent clusters, irrespectively of the number and sort of 
atoms in the cluster. It is only for large rare gas clusters, 
with N ~ 10 3 or more, that we expect a transition to a 
nanoplasma behavior, as it has been found in hydrody- 
namical simulations of such systems [|[ [l2|]. Where and 
how this transition happens will be the subject of further 
studies as well as the connection with enhanced energy 
absorption recently reported for small metal clusters [n3| . 
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